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Abstract: Diffractive excitation is usually described by the Good- Walker formalism for 
low masses, and by the triple-Regge formalism for high masses. In the Good-Walker 
formalism the cross section is determined by the fluctuations in the interaction. In this 
paper we show that by taking the fluctuations in the BFKL ladder into account, it is 
possible to describe both low and high mass excitation by the Good-Walker mechanism. 
In high energy pp collisions the fluctuations are strongly suppressed by saturation, which 
implies that pomeron exchange does not factorise between DIS and pp collisions. The 
Dipole Cascade Model reproduces the expected triple-Regge form for the bare pomeron, 
and the triple-pomeron coupling is estimated. 
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1. Introduction 

Diffractive excitation represents large fractions of the cross sections in pp collisions or 
DIS. In most analyses of collisions low mass excitation is described by the Good-Walker 
formalism [1], while high mass excitation is described by a triple- Regge formula [2,3]. In the 
Good-Walker formalism the state of the incoming projectile is written as a superposition of 
eigenstates to the T-matrix, and the cross section for diffractive excitation is given by the 
fluctuations in the eigenvalues. In the triple-Regge formulation it is instead determined by 
the reggeon couplings to the projectile and the target, and a set of triple-reggeon couplings, 
determined by fits to data (for recent analyses see e.g. refs. [4,5]). The fiuctuations in the 
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pomeron ladder arc here not included in the Good-Walkcr formalism, which therefore 
limits the application to low masses. It is, however, well known that the fluctuations in 
the evolution of a BFKL pomeron are very large [6]. As we will discuss in the following, 
by including these fluctuations it is possible to describe both low and high mass diffraction 
in a uniform way, within the Good-Walker formalism. 

In central pp collisions the interaction is approaching the black limit at increasing 
energy, and therefore saturation effects are very important. The triple-Regge formula would 
violate unitarity and predict a diffractive cross section exceeding the total cross section if 
saturation and multiple pomeron interactions are not included. These are accounted for 
in terms of gap-survival form factors and "enhanced diagrams" , as in the references cited 
above, or as saturation effects in the pomeron flux [7]. In the Good-Walker approach these 
effects are taken into account by reduced fluctuations, when the interaction approaches the 
black limit. 

Data from HERA show a very large cross section for diffractive excitation of the virtual 
photon. In DIS the photon couples initially to a virtual qq pair. To improve the description 
of diffractive excitation of the photon, gluon radiation has been included. The incoming 
virtual photon has been treated as a mixture of qq and qqg states, and the data has been 
fitted to diffractive proton structure functions or parton distributions in the pomeron [8-10]. 
In a description in transverse coordinate space also effects of saturation for small have 
been taken into account [11]. Although important for very small x and small Q^, saturation 
is much less essential in DIS than in pp collisions, which can explain the lack of factorisation 
in the comparison of DIS and pp collisions [12] (see e.g. ref. [13]). 

The eikonal approximation, formulated in impact parameter space, is a formalism 
which efficiently accounts for saturation effects and unitarity constraints in high energy 
reactions. If the colliding particles have a substructure, the eikonal formalism can also 
describe diffractive excitation within the Good-Walker formalism. Miettinen and Pumplin 
[14] suggested that the scattering eigenstates correspond to parton showers, which interact 
via parton-parton scattering. (They also suggested that the partons might be identical to 
quarks and gluons, which at the time were still hypothetical.) The model predicted that 
diffractive excitation is dominantly peripheral, with a maximum for impact parameter 
b « 0.5 fm at = 53 GeV. 

Mueller and coworkers have developed a dipole cascade model in transverse coordinate 
space, which at the same time reproduces leading log BFKL evolution and satisfies s- 
channel unitarity [15-17]. The evolution of the cascade gives dipole chains, which interact 
via gluon exchange. Multiple interactions then correspond to the exchange of multiple 
pomerons. It was pointed out by Mueller and Salam [6] that the dipole evolution contains 
very large fluctuations. This caused a technical problem for their MC simulations, but, as 
discussed below, including the fluctuations in the pomeron ladder gives the possibility to 
treat also higher mass excitations in the Good-Walker formalism. 

In a series of papers [18-21] a generalisation of Mueller's model is presented, which 
includes the following improvements: 

- NLL BFKL effects 

- Nonlinear effects within the evolution 
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- Confinement effects 

- A simple model for the proton wavefunction 

Tliis model describes successfully total and (quasi) elastic cross sections for DIS and pp 
collisions. While taking into account not only fluctuations in the projectile wave function, 
but in the whole evolution between the projectile and the target, the model is also able to 
describe diffractive excitation, not only to low, but also to high masses [20]. Studying a 
collision in a frame, where the projectile is evolved a distance Yp in rapidity, and the target 
a distance It = ^ — it is possible to calculate diffractive scattering where the rapidity 
range of the excited projectile, approximately given by InM-^, is smaller than Yp. (See 
sec. 3.3 for details.) Varying Yp then gives the mass distribution da/dln ~ da/dYp. 

In a similar way it is possible to calculate double diffractive excitation for < 
exp(l^) and M|-^. < exp(y — Yt), where the projectile and target are excited to Mxp and 
Mxt respectively. We note that final states, where the two excited states overlap in rapidity, 
cannot be calculated in this way; in this formalism they are instead included in the inelastic 
cross section. (We want to return to this problem in a future publication.) 

The aim of this paper is to study the nature of the fluctuations in the evolution of 
parton cascades in more detail, in order to understand the relation between the Good- 
Walker and the triple-Regge formalism for diffractive excitation. We will see that within 
the dipole cascade model the Good- Walker mechanism indeed reproduces the expected bare 
pomeron trajectory and the triple-Regge result for diffraction. We will also investigate the 
effects of saturation in more detail, and how the absorptive effects and enhanced diagrams 
correspond to saturation effects in dipole cascade evolutions, and how this describes the 
breaking of factorisation between DIS and pp scattering. 

In the present paper we will not discuss the properties of exclusive final states in 
diffraction, or events with multiple rapidity gaps. We hope to return to these questions in 
future publications. We are also here not discussing the nature of hard diffraction, which 
has been analysed in terms of a hard parton scattering supplemented with extra gluon 
exchange neutralising the colour exchange, together with Sudakov form factors describing 
the gap survival probability (see e.g. ref. [22-25]). 

Section 2 of this paper introduces the Good-Walker formalism and how it can be 
applied to parton cascades. Section 3 summarises the features of the Lund dipole cascade 
model used in our analysis. The nature of the fluctuations and effects of saturation in DIS 
and pp collisions is analysed in section 4, and in section 5 we study the impact parameter 
profile and the t-dependence in pp scattering. The results of the Good-Walker analysis is 
compared to the triple-Regge formalism is section 6, and the bare pomeron couplings are 
estimated. Our conclusions are summarised in section 7. 



2. The eikonal approximation and the Good— Walker formahsm 
2.1 Eikonal approximation 

Diffraction, saturation, and multiple interactions are more easily described in impact pa- 
rameter space. In transverse momentum space the amplitude for two successive interactions 
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is represented by a convolution of the single interaction contributions, which in impact pa- 
rameter space simplifies to a multiplication. 

If the scattering is driven by absorption into a large number of inelastic states n, with 
Born amplitudes \/2/„, the optical theorem gives an elastic Born amplitude 

^ = E/"- (2.1) 

In our notation, where T = 1 — S, these amplitudes are purely real. In the eikonal approx- 
imation multiple interactions exponentiates, and the amplitude 

r = l-e"-^ = 1 -e"^^" (2.2) 

is always satisfying the unitarity constraint T < 1. For a structureless projectile we then 
find: 

datot/dH^ {2T) 
daei/d^b ~ {Tf 

dcTi^ei/dH - (1 - e-S2/„^ ^ datot/cfb - da^^/dH (2.3) 



2.2 Good— Walker formalism 

If the projectile has an internal structure, the mass eigenstates can differ from the eigen- 
states of diffraction. We denote the diffractive eigenstates with eigenvalues T^, and 
the mass eigenstates '^k = J2n ^^kn^n, where the incoming state is given by ^'in = 
The elastic amplitude is then given by (assuming here that cin are real) 

(*i|r|*i) = ^cLTn,= (r) (2.4) 

which implies that 

da,i/d% = (J]cLt„)' = {Tf. (2.5) 
The amplitude for diffractive transition to the mass eigenstate ^'^ becomes 

(2.6) 

n 

which gives a total diffractive cross section (inch elastic scattering) 

da^s/d^b = ^{^,\T\^k){^k\T\^i) = {T^)- (2-7) 
k 

Subtracting the elastic scattering we find the cross section for diffractive excitation 

da^s ,^/d^b = daAss/d% - da^i/dH = (T^) - (r)^ = Vt, (2.8) 
which thus is determined by the fluctuations in the scattering process. 
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Figure 1: (a) An example of a parton (or dipole) cascade evolved in rapidity, (b) The exchange 
of a gluon gives rise to an inelastic interaction, (c) Elastic scattering is obtained from coherent 
scattering of different partons in different cascades, via the exchange of two gluons. (d) Diffractive 
excitation is obtained when the result of the two-gluon exchange does not correspond to the coherent 
initial proton state. Here the dashed lines indicate virtual emissions, which are not present in the 
diffractive final state. 



2.3 What are the diffractive eigenstates? 

As mentioned in the introduction, Miettinen and Pumplin [14] assumed that the diffractive 
eigenstates correspond to parton cascades, which can come on shell through interaction with 
the target. This was also the assumption in our earlier analysis of diffractive excitation 
in [20]. The process is illustrated in fig. ||. Fig. a shows the virtual cascade before the 
collision, and fig. b illustrates an inelastic interaction, where gluon exchange gives a colour 
connection between the projectile and the target. Fig. c shows an elastic interaction, where 
two gluons scatter coherently on the partons in the projectile cascade. It is obtained from 
the projection of the scattered state onto the incoming mixture of different cascades. Fig. 
d, finally, shows the contribution of the scattered state, which is orthogonal to the incoming 
state, and thus corresponds to diffractive excitation. The lines can symbolise gluons in a 
traditional cascade, or dipoles in a dipole cascade. In fig. d the dashed lines corresponds 
to virtual emissions in the cascade, which cannot come on shell via momentum exchange 
from the exchanged gluon pair. 

A similar approach was also used by Hatta et al. [26]. Their analysis was, however, 
limited to relatively low mass excitations. As the authors sought an analytic solution, 
they studied very high energies, where the fluctuations in the pomeron evolution could be 
neglected due to saturation. Thus only fluctuations coming from ordered DGLAP chains 
close to the virtual photon end of the process were included, and therefore it was not 
possible to treat excitation to larger masses. 



3. The dipole cascade model 
3.1 Mueller's dipole model 

Mueller's dipole cascade model [15-17] is a formulation of BFKL evolution in transverse 
coordinate space. Gluon radiation from the colour charge in a parent quark or gluon is 
screened by the accompanying anticharge in the colour dipole. This suppresses emissions at 
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large transverse separation, which corresponds to the suppression of smah k_i_ in BFKL. For 
a dipole {x,y) the probabihty per unit rapidity (Y) for emission of a gluon at transverse 
position z is given by 

dV a 2 (x- yf - 

— — = — a z- -77- -77, with a = . o.l 

This emission imphes that the dipole is split into two dipoles, which (in the large Nc limit) 
emit new gluons independently. The result is a cascade, where the number of dipoles grows 
exponentially with Y . 

In a high energy collision, the dipole cascades in the projectile and the target are 
evolved from their rest frames to the rapidities they will have in the specific Lorentz frame 
chosen for the analysis. The growth in the number of dipoles also implies a strong growth 
for the scattering probability, which, however, is kept below 1 by the possibility to have 
multiple dipole interactions in a single event. The scattering probability between two 
elementary colour dipoles with coordinates {xi,y^) and {xj,yj) in the projectile and the 
target respectively, is given by 2fij, where (in Born approximation) 



fij — f {•^i-iyi\-^j lU j) — Q 



"s f{xi-y^Y{yi-Xj 



log 



{xi - x,jf{y^ -yjf 



(3.2) 



The optical theorem then implies that the elastic amplitude for dipole i scattering off dipole 
j is given by fij. Summing over i and j gives the one-pomeron elastic amplitude 



In the eikonal approximation the unitarised amplitude is given by the exponentiated ex- 
pression 

T{h) = 1 - e-^, (3.4) 
and the total, diffractive, and elastic cross sections are given by the expressions in eqs. 



dp, I, 



3.2 The Lund dipole cascade model 

In refs. [18,19,21] we describe a modification of Mueller's cascade model with the following 
features: 

• It includes essential NLL BFKL effects. 

• It includes non-linear effects in the evolution. 

• It includes effects of confinement. 

The model also includes a simple model for the proton wavefunction, and is imple- 
mented in a Monte Carlo simulation program called DIPSY. Here the NLL effects signif- 
icantly reduce the production of small dipoles, and thereby also the associated numerical 
difficulties with very large dipole multiplicities are avoided. As discussed in the cited refer- 
ences, the model is able to describe a wide range of observables in DIS and pp scattering, 
with very few parameters. 
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3.2.1 NLL effects 

The NLL corrections to BFKL evolution have three major sources [27]: 
The running coupling: 

This is relatively easily included in a MC simulation process. 
Non-singular terms in the splitting function: 

These terms suppress large z-vahics in the individual parton branchings, and prevent 
the daughter from being faster than her recoiling parent. Most of this effect is taken care 
of by including energy-momentum conservation in the evolution. This is effectively taken 
into account by associating a dipole with transverse size r with a transverse momentum 
k± = 1/r, and demanding conservation of the lightcone momentum p+ in every step in the 
evolution. This gives an effective cutoff for small dipoles, which eliminates the numerical 
problems encountered in the MC implementation by Mueller and Salam [6] . 

Projectile-target symmetry: 

This is also called energy scale terms, and is essentially equivalent to the so called 
consistency constraint. This effect is taken into account by conservation of both positive 
and negative lightcone momentum components, p.^ and p-. The treatment of these effects 
includes also effects beyond NLL, in a way similar to the treatment by Salam in ref. [27]. 
Thus the power Aeff, determining the growth for small x, is not negative for large values of 
as- 

3.2.2 Non-linear effects and saturation 

As mentioned above, dipole loops (or equivalently pomeron loops) are not included in 
Mueller's cascade model, if they occur within the evolution, but only if they are cut in 
the Lorentz frame used in the calculations, as a result of multiple scattering in this frame. 
The result is therefore not frame independent. (The situation is similar in the Colour 
Glass Condensate or the JIMWLK equations.) As for dipole scattering the probability for 
such loops is given by a^, and therefore formally colour suppressed compared to dipole 
splitting, which is proportional to a = Ncas/ir. These loops are therefore related to the 
probability that two dipoles have the same colour. Two dipoles with the same colour form 
a quadrupole. Such a field may be better approximated by two dipoles formed by the 
closest colour-anticolour charges. This corresponds to a recoupling of the colour dipole 
chains. We call this process a dipole "swing". The swing gives rise to loops within the 
cascades, and makes the cross section frame independent up to a few percent. We note 
that a similar effect would also be obtained from gluon exchange between the two dipoles. 

In the MC implementation each dipole is assigned one of colours, and dipoles with 
the same colour are allowed to recouple. The weight for the recoupling is assumed to be 
proportional to rfr2/{r^r'l), where ri and r2 are the sizes of the original dipoles and rs 
and r4 are the sizes of the recoupled dipoles. We note that in this formulation the number 
of dipoles is not reduced. The given weight favours the formation of smaller dipoles, and 
the saturation effect is obtained because the smaller dipoles have smaller cross sections. 
Thus in an evolution in momentum space the swing would not correspond to an absorption 
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of gluons below the saturation line k'j_ = Q^(x); it would rather correspond to lifting the 
gluons to higher k± above this line. 

Although this mechanism does not give an explicitly frame independent result, MC 
simulations show that it is a very good approximation. 

3.2.3 Confinement eff"ects 

Confinement effects are included via an effective gluon mass, which gives an exponen- 
tial suppression for very large dipoles [20]. This prevents the proton to grow too fast in 
transverse size, and is also essential to satisfy Froisart's bound at high energies [28]. 

3.2.4 Initial dipole configurations 

Photon wavefunction 

An initial photon is split into a qq pair, and for larger the wavefunction for a virtual 
photon can be determined perturbatively. The well known result has the following form: 



TT 



hh 



^]UQ,r,z) = y^^^^^ej (3.5) 



le 



vr 



z5h+h- - (1 - z)5h-5i+) efKi{r£f) + 5fe+5^+m/i^o(?-e/)} 



with 

£/ = ^z{l- z)Q'^ + m). (3.6) 

Here r is the transverse size of the dipole, and z is the energy fraction carried by the quark, 
A = and + denote the longitudinal and transverse wavefunctions respectively, / denotes 
the quark flavour, and Kq and Ki are modified Bessel functions, e/ is the electric charge 
of the quark in units of the proton charge and rrif the effective mass of the quark. For 
smaller a hadronic component has to be added, as described in more detail in ref. [21]. 
Proton wavefunction 

The internal structure of the proton is governed by soft QCD, and is not possible 
to calculate perturbatively. In the our model it is represented by an equilateral triangle 
formed by three dipoles, and with a radius of 3 GeV~^ « 0.6 fm. The model should be used 
at low X, and when the system is evolved over a large rapidity range the observable results 
depend only weakly on the exact configuration of the dipoles, or whether the charges are 
treated as (anti)quarks or gluons. 

3.3 Application to diffraction 



We now want to apply the Good-Walker result in eq. (|2.8D to the situation where two 
different cascades collide. The elastic scattering amplitude is obtained when T is averaged 
over both the projectile and the target states, while the total diffractive cross section is 
obtained by averaging T^. Thus we have 

da,i/dH = {T)l, (3.7) 
da^n/d^h = {T\t (3.8) 
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^target — X^m ^m^t,r 



target 



Figure 2: Single diffractive excitation with no final state particles in the It range. The virtual 
target evolutions are summed on amplitude level, while the real projectile evolutions are summed 
on cross section level. 



Here the indices p and t indicate averaging over the projectile and target evolutions respec- 
tively. If we average the amplitude over possible evolutions of the target system, we get 
the amplitude representing elastic scattering of the target, If we then square, and average 
over projectile states, we get according to eq. (|2.7| ), the cross section for total diffractive 
scattering of the projectile, while the target is only scattered elastically. Subtracting the 
cross section for elastic scattering of both the projectile and the target gives the cross 
section for single diffractive excitation of the projectile: 

da^roj diScJdH = {{T)X - {T)lt (3.9) 

The process is illustrated in fig. ^. If the expression is calculated in a Lorentz frame in 
which the projectile is evolved a rapidity range Yp, the partons in the projectile cascade are 



confined to the rapidity range y <Yp. The result in eq. (|3.9D includes all cascades limited 
to this range, also those which have no partons close to Yp. This corresponds to all possible 
excitation masses < exp(yp) • IGeV^. By varying Yp it is then possible to calculate the 
differential cross section da^is cx/dM^. Final states with > exp(yp) • IGeV^ are thus 
not included in the cross section in eq. ( |3.9D , in the frame chosen for the calculation. These 
states are in our formalism instead included in the inelastic cross section, because in such a 
frame there is colour exchange connecting the forward- and backward- moving systems. To 
get the full cross section for single diffractive excitation of the projectile, we must do the 
calculation in the target rest frame, where Yp = Y = lns. In the same way it is possible to 
calculate single excitation of the target, by replacing the role of projectile and target. 

The cross section for diffractive scattering of both the projectile and the target is 
obtained by {T'^)pt. This expression includes both elastic scattering and single diffractive 
excitation of the projectile or the target. Subtracting these contributions using eqs. ( ^?^ ) 



and ( |3.9| ), we get the cross section for double diffractive excitation given by 

dann/d'b = {T\, - {{T)X - ((T)^), + (r)^,. (3.10) 
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Figure 3: The fraction of elastic, single diffractive and double diffractive events at 1800 GeV as 
function of interaction frame. For single diffractive excitation the figure shows masses of the excited 
projectile integrated over M|p < exp(yp) • IGeV^. For double diffraction the projectile and target 
masses are integrated over M|p < exp(rp) • IGeV^ and M^^ < exp(yt) • IGeV^ respectively, with 
Yp + Yt = Y = ln(1800^) « 15. (For details see the main text.) The two lower error bands are single 
diffractive excitation and elastic cross section estimated from CDF data [29,30]. The top area is 
the sum of the two, thus not including double diffraction. 

This expression gives the cross section for < exp(yp) • IGeV^ and M^^ < exp(yt) • 

IGeV^, where Yp+Yt equals the total rapidity range Y. As was the case for single diffractive 
excitation, events with excitation to larger masses are in this formalism included in the 
inelastic cross section. For single diffraction it was possible to include excitation of e.g. 
the projectile to all masses by performing the calculation in the target rest frame. This 
is not the case for double diffraction. Even if we change Lorentz frame, we can never 
include events where the two excited states overlap in rapidity. Those states will always be 
included in the inelastic cross section. (Thus although the total and elastic cross sections 
have to be independent of the Lorentz frame used, only the sum of the cross sections for 
inelastic scattering and diffractive excitation is frame independent.) 

The results from MC simulations of single and double diffractive excitation were pre- 
sented in ref. [20], in good agreement with data from HERA and the Tevatron. Fig. ^ shows 
the diffractive cross sections for pp collisions at 1800 GeV. In this figure the projectile is 
evolved over Yp units of rapidity, and the target over Yt = Y — Yp units, setting the limits 
for the diffracted masses to < exp{Yp) ■ IGeV^ and, in case of double diffraction. 



We will in the next two sections study how the results follow from the nature of the 
fluctuations causing the excitations, and how the fluctuations are suppressed by saturation 



< exp(Y - Yp) ■ lGeV2. 
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Figure 4: Probability distribution, P{F), for the one-pomeron araplitude F in DIS, represented 
by a dipole with size r = 1/Q, for = 14 QeV^ and W = 220 GeV (left) and W = 1000 GeV 
(right), b is in units of GeV^^. The dotted lines are fits of the form in eq. (4.1). 



effects. We also note that in this approach the effective triple-pomeron coupling is fixed 
by the constraint, that it is the same dynamics that determines both the coupling between 
the three pomeron ladders in fig. |2|, and the evolution within the individual ladders. The 
relation to the triple-Regge formalism will be discussed in sec. 6. 

4. The nature of the fluctuations and eff'ects of saturation 
4.1 7*p scattering 

The photon wavefunction in eq. ( |3.6| ) is divergent for small dipole sizes, which means that 
infinitely many small dipoles are created with infinitely small cross sections. To illustrate 
the fluctuations in the dipole cascade we show in fig. ^ MC results for the probability 
distribution, P{F), for the one pomeron amplitude F in eq. ( ^.iD for a dipole with a fixed 
size r = l/Q at a fixed impact parameter b. The distribution P{F) is here defined so that 
P{F)dF is the probability for the formation of a pair of a projectile and a target cascade, 
for which the Born amplitude F = ^2 fij li^s between F and F + dF. The calculations are 
performed in the hadronic cms, which implies that the diffractive masses are integrated 
over the range Mjr < \^W^ ■ 1 GeV. 

As seen in fig. ^, the probability distributions can for all 5- values be well approximated 
by a power spectrum 

P{F)^AF-P, (4.1) 

with a cutoff for small F-values. These approximations are shown by the dotted lines. The 
two parameters A and p are tuned to fit the MC results for different values of the energy 
W, dipole size 1/Q, and impact parameter b. (The cutoff is then adjusted to satisfy the 
normalization condition J P{F)dF = 1.) As we will see below, it is particularly interesting 
to note, that the fitted value for the power p is independent of the impact parameter b. It 
varies, however, slowly with and W as can be seen in table |l]. 

The cross sections obtained from these distributions can most easily be estimated from 
the approximation in eq. ( |4.1| ). We see in fig. ^ that the Born amplitudes are generally 
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Table 1: The values of the parameters in the fits to P{F) for ^*p, eq. (4.1), and pp, eq. (4.4), for 
some sample energies and 6-values. The power p is independent of impact parameter in 7*p, while 
for pp it is the exponential suppression a that does not depend on b. 



small, which implies that miitarity effects are small, and T = 1 — ~ F. We also note 
that the widths of the distributions are large, which means that (T)^ can be neglected 
compared to (T^). The approximation in eq. ( [4.1[) then gives the result 



dotot ^ = 2A I (1 - e-^)F-PdF = -2^r(l -p); 
Jo 

-Vt = (T^) - {Tf « (T^) = (4.2) 



Aj{l- e'^F-^dF = (1 - ^) X 2{T). 



From these results we note that the ratio dfJdifr cx/^^ctot = ^/(2(r)) depends only on the 
value of the parameter p. As we have found that p is independent of the impact parameter 
for fixed W and Q^, we can integrate over b, and find 

atot 2{T) ^ 

Thus the parametrisation in eq. ( ^ ) gives crdifi ex/ftot ~ 0.18 for = 14GeV^ falling to 
~ 0.13 at = 50GeV^. Although the simple parametrisation overestimates the result of 
the MC, it gives a qualitatively correct result. 

For a virtual photon in DIS the fluctuations will be further enhanced by adding the 
fluctuations in the photon wave function, but this will not alter the conclusions presented 
above. 

4.2 pp scattering 

The corresponding Born amplitude distributions in pp collisions are shown in fig. ^ for 
W = 100 and W = 2000 GeV and different 6-values. We note that here the interaction 
probability is large, which implies large saturation effects. The distributions can be well 
approximated by Gamma functions of the form 

P{F) = AFPe-"^. (4.4) 

The distributions have two parameters, p and a, which are tuned to the MC results for 
different values of energy and impact parameter. The parameter A is then fixed by the 
normalisation condition. The result of the fit is shown in table |l], and we note here that 
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a is essentially independent of the impact parameter, but falling with energy. For a fixed 
energy, the decrease in (F) for more peripheral collisions is related to a decrease in p for 
larger 6- values, pushing the distribution to smaller values of F. For fix b the parameter p 
also grows with increasing energy, reflecting the larger interaction probability. 
The parametrisation in eq. (|4.4|) gives 



a 

Vf 1 



2(F) 2a 



0.35 for W = 100 GeV, (4.5) 



where Vp = {F'^) - (F)^ is the variance of F. Thus we find also here that the ratio between 
the variance and the average of the Born amplitude is independent of b. We note that 
this ratio is large, and similar to the result for 7*1? collisions at lower Q^-values. Thus, 
without saturation we would have a correspondingly large value for da^is ex/dctot also in 
pp collisions. 

However, as the one pomeron amplitude (F) is large in pp scattering, unitarity cor- 
rections are very important. The probability distribution for the unitarised amplitude, 
F(r;6), with T = 1 - e"^, is shown in fig. | for = 100 and 2000 GeV and different 
6- values. We see that for the central collisions the distributions are very peaked close to the 
unitary limit T = 1. This reduces the fluctuations very strongly. For the parametrisation 
in eq. ([4. 4]) the average and the variance for the distribution in T are also easily calculated, 
and given by 



(T) = 1 



'a + l' 



Vt = i^r' - i^f'^'- (4.6) 
a + 2 a + 1 

We see that at high energies and central collisions, where the Born amplitude (F), and 
thus also the parameter p, become large, (T) will approach 1 and Vr will go towards 0. For 
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Figure 6: Probability distribution, P{T), for the full amplitude T for pp collisions a.t W = 100 
GeV (left) and W = 2000 GeV (right), b is in units of GeV"^ The dotted lines correspond to fits 
of the form in eq. (4.4). 



central collision at W = 100 GeV the ratio of diffractive events Vt/(2(T)) is about 0.035, 
a factor 10 lower than without unitarisation. Therefore in central collisions diffractive 
excitation is suppressed, and diffractive scattering is dominantly elastic. 

The large effect of saturation in pp collisions has also the effect that factorization is 
broken when comparing diffractive excitation in DIS and pp collisions [12]. 



5. Impact parameter profile and t-dependence in pp-collisions 

In the previous section we showed the amplitude fluctuations for different impact parameter 
values. We will here study the 6-dependence, and the corresponding t-dependence, in more 
detail. As mentioned above, diffractive excitation is small in central pp collisions as (T) 
is approaching 1. In highly peripheral collisions both (T) and Vr are small, again giving 
little diffractive excitation. Therefore diffractive excitation is dominated by moderately 
peripheral collisions, where (T) ~ 0.5 and {F) ^ 1. The 6-dependence of (Ttot/2, o"ci, and 
Cdiff ex in the MC is shown in fig. 0. 

As pointed out also in earlier analyses {e.g. in refs. [14,31]), this implies that diffractive 
excitation in pp collisions appears in a ring with a radius which grows slowly with energy. 
In a purely perturbative calculation with massless gluons, the total cross section will grow 
very fast due to the formation of very large dipoles, and eventually violate Froisart's bound. 
However, as demonstrated by Avsar [28] , the inclusion of confinement effects via a massive 
gluon (as in the simulations described above) implies that very large dipoles are suppressed, 
and the black disk radius grows proportional to Ins. This means that the total and elastic 
cross sections grow like In^ s for very large energies. In addition the results in [28] show that 
the slope, when the interaction drops from black in central to white for more peripheral 
collisions, is approximately constant with energy. Thus the width of the ring with large 
diffractive excitation is approximately constant at high energies. Consequently the cross 
section for diffractive excitation will for very large energies grow proportional to the radius 
of the ring, i.e. proportional to Ins. 
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Figure 7: Impact parameter distributions from tlie MC for (T) — {da tot I <^V) j 2, (T)^ — daci/d^b, 
and Vt = dadiSc^/d'^b in pp collisions at W ^ 100, 2000, and 14000 GeV. b is in units of GeV"^ 



In our model the interaction is driven by absorption into inelastic channels, and with 
our definition, where S = 1 — T, the imaginary part of the amplitude T is neglected. It 
is therefore straight forward to take the Fourier transform and calculate the i-dependence. 
The result for the differential elastic pp cross section was presented in ref. [21], and the 
result for single diffractive excitation is shown in fig. |^. Figure |^ shows the result at 546 
GeV, together with an extrapolation of a fit to UA8 data [32], normalised to the model 
result. The UA8 data cover only the range 0.8 < |t| < 2.5 GeV^, but we note that the 
t-slope in the model agrees well with this fit. For comparison also the elastic cross section 
is included in the figure, and we see that the diffractive slope is significantly smaller than 
the slope in elastic scattering. This is a consequence of the larger 6-values for diffractive 
excitation shown in fig. 0. 

Figure ^ shows how the differential cross section for diffractive excitation is varying 
with energy. We see that the increase is quite slow, as a result of saturation and unitarity 
constraints. Note that for low t, the energy dependence is stronger, due to the growth of 
the radius of the ring, while the energy dependence for high t is much slower, showing that 
the width of the ring is almost energy independent, in agreement with ref. [28]. 

To further illustrate the effect of saturation we show in fig. |8^ also the t-dependence of 
the single diffractive cross section obtained from the unsaturated Born amplitude. We see 
that saturation reduces the cross section by roughly a factor 25 at 2000 GeV. The slope is 
less affected, but the suppression for small 6-values implies that the t-dependence deviates 
more from a pure exponential, when saturation is included. 

6. Relation Good- Walker — Triple-Regge 

In this section we will discuss the relation between the results using the Good-Walker 
formalism described above, and the triple- Regge formalism. In this comparison we want 
to study the contribution from the bare pomeron, meaning the one-pomeron amplitude 
without contributions from saturation, enhanced diagrams or gap survival form factors. 



-15- 




Figure 8: (a) t-dependence for the single diffractive cross section at 546 GeV from the MC, together 
with a fit to UA8 data. The elastic cross section is included for comparison, (b) MC results for the 
energy dependence of the ^-distribution for single diffractive excitation. The effect of saturation is 
demonstrated by the dotted line, which shows the result at 2000 GeV without saturation, scaled 
by a factor 0.1. 



We want to see if the fluctuations in the dipole cascades reproduce the powerlike energy 
dependence expected in the Regge formalism. 

When s, M|, and s/M| are not small, pomeron exchange should dominate. If the 
pomeron is a simple pole we expect the following expressions for the pp total and diffractive 
cross sections: 



dt IQir' 



rl ^ / \ 2(aW-l) 

^dtd{Ml) 16tt'^ ^ ^^^^ ^ \MlJ ^ ^ ^ 

Here a{t) = 1 + e + a't is the pomeron trajectory, and and g3p{t) are the proton- 
pomeron and triple-pomeron couplings respectively. (We have here omitted the scale sq in 
the powers (s/so)° or (M|/so)". This scale is in the following is assumed to be 1 GeV .) 

The results of the MC for the total, elastic, and single diffractive cross sections are 
shown by the crosses in fig. |^. The elastic and diffractive cross sections are integrated 
over t and M|. The single diffractive cross section is calculated in the total cms, which 
corresponds to an integration over masses in the range < GeV, and it corresponds 

to excitation of one side only. We see that the result indeed has the powerlike increase 
with energy, which is characteristic for a Regge pole. We also note that in the one-pomeron 
approximation the elastic cross section is larger than the total for ^/s > 15 GeV. 

If we assume a simple exponential form for the proton-pomeron coupling, /3(t) = 
/3(0) exp(6o,ei '^^^ integrate the elastic cross section in eq. ( |6?l] ) over t, and obtain 

2 

^ei = TTr^TT' ^i*"^ ^(*) = ^o,ei + 2q' In s. (6.2) 
167rB{s) 
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Figure 9: The total, elastic and single diffractive cross sections in the one-pomeron approximation. 
The crosses are from the dipole cascade model without saturation, and the lines are from a tuned 
triple Regge parametrisation. 



Besides the shrinking of the elastic peak, the pomeron slope a' also gives a logarithmic 
correction to the powerlike increase of the elastic cross section. A consistent fit to both 
quantities is obtained for a' = 0.2 GeV~^. In fig. |^ the lines are obtained from the 
expressions in eqs. |6.2D with the parameter values 



a(0) = l + e= 1.21, a' = 0.2GeV"^ 

= ^2(0) = 12.6mb, 6o,ei = 8GeV-^ g3p{t) = const. = 0.3GeV" 



(6.3) 



We have here assumed a constant triple-pomeron coupling, and we see that the MC results 
in fig. 1^ are very well reproduced by this fit. 

The t-dependence of elastic scattering and diffractive excitation, shown in fig. are 
however not pure exponentials, as assumed in the fit above. (For diffractive excitation a 
minor deviation from a pure exponent originates from the integration over M^.) We want 
to study this dependence in some more detail, and are here in particular interested in the 
t-dependence of the triple-pomeron coupling g3p{t). A very close fit to dcTd/dt is obtained 
for (with t measured in GeV^) 



/3^t) = /34(0)exp 



lot 



1 - l.St 



(6.4) 



which cannot be distinguished from the MC result in fig. 10. Inserting this fit into the 



expression for the diffractive cross section integrated over and assuming a constant 
triple-pomeron coupling equal to 0.27 GeV~^, gives the thin dotted line in fig. 10. We 
see that this is quite a good fit. It is slightly less steep for |t|-values below 0.1 GeV^, but 
deviating less than 10% from the result of the model. We also note that these modifications 
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Figure 10: i-dependence for elastic and single diffractive excitation without saturation effects, at 



1800 GeV. A very close fit to da^i/dt is obtained from the proton-pomeron coupling in eq. (6.4). 
Including a constant triple-pomeron coupling in the expression for diffractive excitation gives the 
thin dotted line, which gives a good description of the model result. For comparison also a pure 
exponential fit to the elastic cross section is included. 



of the t-dependence for the elastic and diffractive cross sections do not modify the good 
fits to the integrated cross sections in fig. ^ 

We here want to make the following comments: 

• Comparison with perturbative QCD 

The expressions in eq. (6.1) correspond to a pomeron which is a simple pole. This 



is not the case in perturbative QCD. In the LL approximation the pomeron is a cut 
in the angular momentum plane, which gives logarithmic corrections to the proton- 
pomeron coupling: /3(0) ~ l/(lns)^/'^ while /3(t) ~ l/(lns)^/^ for t 7^ [16,33]. The 
different s-dependence when t is equal to, or different from, zero is associated with a 
cusp in the t-dependence at t = 0. However, at these small i-values perturbative QCD 
is not applicable, and non-perturbative effects are important. It was also pointed out 
by Lipatov [34], that a running coupling can modify the cut to a series of poles. (Note 
that a running coupling is included in our model simulations.) A strong increase in 
the elastic cross section at very small t-values is not seen in the experimental data, 
and also not present in our result in fig. |l^. An extra factor of l/(lns)^/^ in /3(0) 
and I3{t) would give an equally good fit to the results in fig. ^, provided the pomeron 
intercept is increased to a(0) = 1.25. It is, however, not possible to find a good fit if 
f3{t) is proportional to l/(lns)^/^, inside a f-region essential for the integrated elastic 
cross section. 

In LL perturbative QCD also the triple-pomeron coupling has a singular behaviour 
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at t = [16,33]: 

1 1 

(lnM2)i/4(lns/Mi)3/4 7=! ^^'^^ 

We saw above that our results were well reproduced by a constant triple-pomeron cou- 
pling, only slightly underestimating the slope for small t- values below \t\ = 0.1, and 
our analysis does not support a triple-pomeron coupling with a strong t-dependence. 
A slowly varying triple-pomeron coupling is also in agreement with early analy- 
ses [35,36]. We note, however, that the magnitude of the triple-pomeron coupling 
agrees (within the large uncertainties) with the perturbative estimate in ref. [33], 
which in our notation corresponds to irg^p ~ 0.2 — 1.7GeV^^ for 0.25 < \t\ < 4GeV^. 

• Comparison with other analyses 

We should also note that the bare pomeron is not an observable. Here we have 
included NLL effects and confinement in the evolution, but not nonlinear effects from 
saturation or multiple collisions. The effect of nonlinearity will appear differently if 
one first adds NLL effects and then compares the results with and without saturation, 
as compared with an approach where the LL result is compared with and without 
saturation, before NLL effects are included. Therefore the bare pomeron may look 
different, depending upon the scheme used to remove the nonlinear effects. 

Keeping this in mind, we want to compare our result in eq. (^^) with some recent 
more traditional Regge analyses. As examples the bare pomeron in the analysis by 
Ryskin et al. [4] has three components with different dependence on the impact pa- 
rameter, or t, in order to mimic the branch cut structure of the pomeron singularity. 
The three poles have the same intercept equal to a(0) = 1.3, and quite small slopes. 
The dominant component has a' = 0.05 GeV~^, and the slope of the other compo- 
nents are even smaller. Ostapchenko [37] finds in an approach with two pomerons 
a(0) 1.35, a' ^ O.OSGeV-^ and a(0) ^ 1.15, a' 0.14GeV^2 for the hard and 
soft pomerons respectively. Gotsman et al. [24] find in an analysis with a single 
pomeron a(0) = 1.335 and a' = O.OlGeV"^. In another fit with a single pomeron 
pole Kaidalov et al. [38] find a(0) = 1.12 and a' = 0.22 GeV~^. We can also compare 
with the results by Goulianos [7,13], who in a formalism with a renormalised pomeron 
flux flnds the values a(0) ~ 1.11 and a' = 0.26 GeV~^. We see here that the way 
saturation is taken into account can have a large effect on the result for the bare 
pomeron. We also note that our result lies somewhere in between these examples. 

In conclusion we see that the Born amplitude in our dipole cascade model indeed 
reproduces the triple-Regge formula for a bare pomeron pole with a(0) = 1.21, a' = 0.2, 
and an almost constant triple-pomeron coupling. The i-dependence of both the elastic and 
the diffractive cross sections is close to an exponential. There is no indication for more 
dramatic variations, like those obtained in LL perturbative QCD, where the pomeron is a 
cut singularity in the angular momentum plane. Here one also expects /3(0) ~ l/(lns)^/^ 
while j3(t) ~ l/(lns)3/^ for t 7^ 0. Our result would be consistent with a proton-pomeron 
coupling proportional to l/(lns)^/^, if the intercept is increased to 1.25, but not with 

~ l/(lns)3/4. 
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7. Conclusions 



Diffr active excitation represents a large fraction of the cross section in pp collisions or DIS. 
In the Good-Walker formalism diffractive excitation is determined by the fluctuations in the 
scattering amplitude. In traditional applications this formalism has been limited to small 
mass excitation, while excitation to high masses has been described in the triple-Rcggc 
formalism, introducing a set of parameters for the pomeron couplings. It was demonstrated 
by Mueller and Salam [6] that BFKL evolution contains large fluctuations, and in this 
paper we demonstrate that, by including these fluctuations in the analysis, it is possible to 
describe diffractive excitation to both low and high masses. 

The Lund Dipolc Cascade model, implemented in the DIPSY MC, describes success- 
fully the total, elastic, and diffractive cross sections in DIS and pp collisions [18,19,21]. In 
this paper we show how the fluctuations in the BFKL evolution can reproduce diffractive 
excitation to low and high masses within the Good-Walker formalism, with parameters de- 
termined only from the total and elastic cross sections. In DIS at HERA the fluctuations 
give a diffractive cross section of the order of 10%, and saturation has a relatively small 
effect on the result. However, in pp collisions unitarity constraints and saturation reduce 
the fluctuations, when the scattering approaches the black disc limit. Therefore diffractive 
excitation in high energy pp collisions is dominated by peripheral collisions, and pomeron 
exchange in DIS and pp collisions does not factorise. 

In the triple-Regge formalism, saturation effects are included in terms of "enhanced 
diagrams" , gap-survival form factors or saturation effects in the pomeron flux, which reduce 
the effect of the "bare pomeron". In this paper we have also studied the effective bare 
pomeron, corresponding to the result obtained when non-linear effects are not included. 
We see that the result indeed is well described by a bare pomeron pole, with a{0) = 1.21 and 
a' = 0.2 GeV^^. The triple-pomeron coupling is fixed by the couplings in the BFKL ladder, 
and the results are well reproduced by a constant coupling gsp Ri 0.3GeV~^. Although 
our model is based on perturbative QCD, with non-pcrturbativc effects introduced only 
in the incoming proton wavefunction and confinement effects via an effective gluon mass, 
this rcsTilt contrasts to LL BFKL, where the pomeron is a cut singularity and the pomeron 
couplings have strong t-dependencies. 

In this paper we have not discussed the properties of exclusive final states in diffractive 
excitation (apart from the mass distribution). We hope to return to this problem and to 
hard diffraction in future work. 

The main results can be summarised as follows: 

- The fluctuations in a BFKL ladder arc large. Taking these fluctuations into account 
as in the DIPSY MC, it is possible to describe diffractive excitations to both low and high 
masses within the Good-Walker formalism. 

- Saturation effects are small in DIS, but large in pp collisions. Therefore diffractive 
excitation is a peripheral process in pp scattering, and factorisation of pomeron exchange 
is broken. 

- The result of the Good-Walker formalism reproduces the triple-Regge result for 
diffractive excitation, with a bare pomeron pole with a(0) = 1.21, a' = 0.2GeV~^, and an 
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almost constant triple-pomeron coupling g^p 0.3 GeV ^. 
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